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0\ . Abstract 

We consider dissipative relativistic fluid theories on a fixed flat, globally 
,^ ■ hyperbolic, Lorentzian manifold (R x T^,gab). We prove that for all initial 

data in a small enough neighborhood of the equilibrium states (in an ap- 
propriate Sobolev norm), the solutions evolve smoothly in time forever and 
decay exponentially to some, in general undetermined, equilibrium state. To 
prove this, three conditions are imposed on these theories. The first condition 
K^ I requires the system of equations to be symmetric hyperbolic, a fundamental 

Q ■ requisite to have a well posed and physically consistent initial value formula- 

(^ . tion. The second condition is a generic consequence of the entropy law, and 

CS I is imposed on the non principal part of the equations. The third condition is 

r^^ ' imposed on the principal part of the equations and it implies that the dissipa- 

0^ ■ tion affects all the fields of the theory. With these requirements we prove that 

Q i all the eigenvalues of the symbol associated to the system of equations of the 

^^1 fluid theory have strictly negative real parts, which in fact, is an alternative 

^ ■ characterization for the theory to be totally dissipative. Once this result has 

been obtained, a straight forward application of a general stability theorem 






% 



due to Kreiss, Ortiz, and Reula, implies the results above mentioned. 
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I. INTRODUCTION 

In recent years there has been a substantial improvement on our understanding on 
how a proper description of dissipative fluids can be incorporated in the framework 
of the theory of relativity. Dissipative relativistic fluid theories satisfying an en- 
tropy law, and having a well posed (symmetric hyperbolic), and causal initial value 
formulation have been presented [^ ^ ^. 

An important result on the physical meaning of all these hyperbolic theories was 
obtained recently 0, §]• It was shown that certain constitutive relations between 
the variables in the hyperbolic fluid theories, which have a clear physical meaning, 
approach in their time evolution the values predicted by the simplest covariant 
generalizations of the Navier Stokes fields. This result is based in a fundamental 
hypothesis namely, that the solution of the hyperbolic fluid field equations exists 
and remains smooth and small during a long enough time interval, such that the 
relaxation to near Navier-Stokes behavior occurs. This is a very important check 
for these fluid theories, since at microscopic scales they are substantially different 
from the usual Navier-Stokes theories in the following sense: for the last ones, one 
expects to have smooth global solutions for all smooth initial data sets, as has been 
proved in lower dimensions , while for the former ones, one expects the developing 
of discontinuities in the form of shock waves for crispy enough initial data. Thus, we 
can only hope to find a neighborhood of equilibrium data for which global solutions 
exist and so, where the departures from Navier-Stokes are uniformly small. 

The purpose of the present work is to look for conditions under which the funda- 
mental hypothesis mentioned above is satisfied. To this end, we apply a theorem 0] 
which is a generalization to the case of partial differential equations of the Ljapunov 
stability theorem for ordinary differential equations. This generalization holds for 
hyperbolic systems such that the eigenvalues of their associated symbols have all 
strictly negative real parts which, as we shall see, is the case for the hyperbolic 
dissipative fluids. 

In order to apply this general stability theorem, three conditions are imposed on 
the dissipative hyperbolic fluid theories. The first condition requires the symmetric 
hyperbolicity of the system of equations; a fundamental requisite to have a well posed 
and physically consistent initial value formulation. As we shall see, the symmetry 
is an automatic property of these theories, but the hyperbolicity has to be required. 
We also include in this condition that the spacetime manifold is R x T^ with a 
flat metric (which is the more restrictive assumption). The other two conditions 
are of a generic type in the sense that all fluid systems, except for very specific 
and isolated ones, satisfy them. Specifically, the second condition requires that the 
nonprincipal part of the system of equations, which is responsible for the dissipation, 
satisfy certain negative-semidefiniteness condition. This condition assures that all 
perturbations to an equilibrium state, which are not tangent to the equilibrium 



submanifold, do dissipate towards equilibrium, and this is manifested by the fact 
that they make a positive definite contribution to the entropy. The third condition 
is a requirement on the principal part of the system of equations, and it means that 
the presence of dissipation affects all the fields of the theory, in the sense of not 
allowing for a decoupled set of fields with its own evolution not being driven by 
dissipation. Both, the second and the third conditions, have already been required 
in the literature with the aim of characterizing the equilibrium states. 

These conditions allow us to apply the theorem proved in |^, which implies not 
only the global existence of solutions, but also their exponential decay to equilibrium 
for initial data near enough, in some appropriate norm, to equilibrium data. 

The plan of the paper is as follows: In Sec. II we briefly introduce the fundamen- 
tal aspects of the fluid theories, and state in detail the character of the conditions 
we impose on them. In Sec. Ill we state and prove our main result. Finally, in Sec. 
IV we present the conclusions. 

II. DISSIPATIVE RELATIVISTIC FLUID THEORIES 

In this section, following [Q, we introduce the dissipative relativistic fluid theories. 
After introducing them, we describe the properties of these fluids needed to prove 
stability. 

We assume that a fluid state is characterized by a finite collection of space-time 
tensor fields. Let i^^ denote these fields, where upper case indices stand for the 
entire set of tensor indices represented in this collection of fields. So we refer to 
ip^ as a point in the space of fluid states S. Lower case indices will denote space- 
time indices. Repeated indices indicate contraction as usual in the abstract index 
notation. We restrict consideration to the fluid theories in which the field equations 
take the form 

MlB^aV"" = -Iabv''- (1) 

where M^^ and Iab are smooth functions of the fluid fields ip^ and the space-time 
metric. 

We say that system (|l]) is symmetric hyperbolic if M^^ = M?j^^^ and there exists 
a timelike-future directed m" such that Nab = ~UaM'^^ is positive definite. This is 
a sufficient condition to have a well posed initial value formulation. 

On physical grounds, to have as maximum speed of propagation the speed of 
light, the stronger condition of causality, that is, —UaM^^ is positive definite for 
all future-directed timelike vector m", is usually required. This means that the 
characteristic surfaces of system (|l|) are inside the null cone given by the space- 
time metric. To prove global existence and decay only the weaker condition of 
hyperbolicity is needed. 

We analize now the structure of the equilibrium states. Following p|, |^ we say 
that a fiuid state ip"^, solution of the dynamical system of equations, is a strict equi- 



lihrium state if its time reverse is also a solution. We denote a strict equilibrium 
state by y?^, and assign a sub index zero to any tensor evaluated at an strict equilib- 
rium state. This definition implies that loABfo = 0. More generally, we say the ip^ 
is a momentary equilibrium, state if Iab^^ = 0- These states are called equilibrium 
states because their entropy production vanishes, and momentary because if this 
condition holds at a certain time, it will not necessarily hold in subsequent times 
(see m). Every fluid state can be written as (p^ = ip^ + rj"^, where Iab^^^ = 0, and 
Tj"^ is such that IabV^ = implies rj^ = 0. We call ip^ the momentary equilibrium 
part of the fluid state. 

Below we state the three conditions imposed on these fluid theories that will be 
used to show global existence and exponential decay. 

1. The fluid system of equations (|lD is symmetric hyperbolic and the space-time 
is (R X T^,gab), where T^ denotes a three dimensional torus and Qab is a flat 
metric. 

2. The tensor Iqab niust be symmetric and positive semidefinite. 

3. The map J^k : 5^ — ^ 5* defined by J-'xi'ip'^) = KaMQAB 'ip^, is injective for 
all space-time vectors K"^ ^ 0; where S* denotes the dual of the space of fluid 
states and S^ the subspace of momentary equilibrium states. 

The first condition is the more restrictive and more work has to be done in order to 
weaken it. It would be interesting to treat physically relevant boundary conditions 
and arrive to similar results. The other two conditions are not very restrictive and 
all fluid systems, except for very specific and isolated ones, satisfy them. 

The second requirement ensures that the effect of Iab in equation (|l|) is to 
dissipate, in the sense offending to move non equilibrium states towards equilibrium 
as time grows. This is just a little stronger than the entropy condition in these fluid 
theories, that requires the entropy source to be non negative. This stronger condition 
was already considered in 0. 

The third requirement is on the principal part of the equations. It implies that 
dissipation affects all the fields of the theory, in the sense of not allowing the ex- 
istence of a decoupled subset of fields with its own evolution not being affected by 
dissipation. This requirement turns out to be equivalent (see appendix A), at least 
for the case of divergence type fluid theories, to one assumed in the literature to 
show that the only possible strict equilibrium states are the constant states ^. 

III. GLOBAL EXISTENCE AND EXPONENTIAL DECAY 

In the theorem below we present our main result about global existence in time and 
decay to strict equilibrium. 



Theorem 1 Consider the Cauchy problem for system ^, corresponding to a hy- 
perbolic divergence type fluid theory satisfying conditions 1-3. If the initial data is 
smooth and close enough in a Hp{T'^) Sobolev norm (p > 5) to the data correspond- 
ing to some strict equilibrium solution, then the solution is smooth, exists globally in 
time and decays exponentially to some strict equilibrium solution in the H^ norm. 

Proof: We define Nab = -UaM%B, N%b = -qtM% where qab = Qab + UaUb is 
the 3- metric in each hypersurface orthogonal to m" (assuming u"'Ua = — 1). Assume 
that the initial data is close to some strict equilibrium state ip^, then the fiuid state 
is ip^ = (Pq + 6ip^. For the variable S(f^, the system (^ becomes 

NABu'^VaS^'' = NXb'^J^'' - IabS^"". (2) 

Here the tensors Nab, ^abi ^^^ Iab are thought as functions of dip"^. The condition 
(1) of Sec. II on the fiuid theories implies that we can choose Cartesian coordinates 
{t, a;-'} on the space-time manifold Rx T^ and a constant u"" such that d/dt = u°'Va- 
Then the system (|^) becomes 

NABp^' = Ni,^ - Iab5v^. (3) 

In order to study solutions near 5(p^ = 0, it is convenient to introduce a parameter 
e to control the smallness of initial data. Thus, dip'^it = 0) = ef^{x^) and the 
solution shall be written as Sip"^ = ev^. As the tensors Nj^^ and Iab are smooth 
functions of ev^, they can be written as 

KB = KAB + eNij,j,, 
Iab = Iqab + ^Iiab- 
With this decomposition the Cauchy problem for (^ is 

Nab^ = (n^ab + eNij^j,) -^ + (loAB + ehAB)v''. (4) 

v^it = 0) = /^(x^). 

As there are periodic boundary conditions on the space coordinates {x-'}, v^ can be 
expanded in Fourier series, 

v^= Y, v^{y,t)e'^-'', 

where Q is the discrete set of Fourier frequencies. 

We want to apply the stability theorem proved in that, for completeness, 
we state in the appendix B. To prove the theorem |l|, we consider the eigenvalues 
problem 

XNoAB^"" = (ikjN^AB - Ioab)^''- (5) 



Then, as explained in the appendix B, we only need to verify the following conditions, 
(i) There is a constant 6 > such that the eigenvalues A(/c-') satisfy Re{A} < —S 

for all y G n, y ^ 0. 

(ii) For y = 0, A(0) < -5 or A(0) = 0. 

(iii) As linear maps, the null space of Iiab contains the null space of Iqab- 
Conditions (ii) and (iii) are satisfied for these fluid theories, since the kernel of Iab 
is of constant dimension (see P, 0), and because condition (2) of Sec. II holds. 

To prove (i), let y be different from zero, then k = Jykj > const. > 0. The 
eigenvalue problem (|^) can be written as 

jNoAB + ^jN'oAb)^"" = {^NoAB + \loAB ' ^jN^Ab)v''- 

Defining Ka = —{\/k)ua + i{kj/k)qi, this can be written as 

KaM^ABi''' = {^NoAB + lloAB ' i^Ni^j,)f^^ . (6) 

Since KaM^j^^ is injective, by condition (3), and a smooth function of kj/k, there 
is a constant c > such that 

D D 

Then, contracting (|^) with itself 

Perturbation theory of linear operators tells us that \{y)/k is uniformly bounded 
(see 0), then we get for some positive c' 

NoBci^ ^^ < c' NoBcff- 
This inequality, together with the positive definiteness of Nqab implies, 

iVoAB^V<5Aro^B^V, 5>0. (7) 

Now, contracting (y) with (p and taking real part 

RelAjA/oAsV^ ^ = -loAB'f V^ = -IqabV i] 

< -~6NoABtv'' 

. 5 -A^B 

< -- NqabV ^ ■ 

c 



Here, 6 > exists because of the negative definiteness of Iqab in the direction of fj^ 
(by condition (2)), and we have used (0) in the last hne. We have thus shown that 

Re{X{k^)} < -5 < 0, with 6 = ->0 and A;^V 0. 

IV. CONCLUSIONS 

In this work we have proved global existence and exponential decay to strict equi- 
librium of solutions, corresponding to initial data in a small enough neighborhood 
of strict equilibrium, for a generic dissipative relativistic fluid theory. 

This result, in particular, verifies a fundamental hj^othesis of previous works 
[^, ^] namely, the existence of solutions during a long enough time interval. The 
closeness of initial data to strict equilibrium is a natural limitation in the sense 
that, for large data, shock waves develop, which are a widely observed phenomena 
in nature. This occurs because these fluid systems are genuinely non-linear. 

There are three questions related to the techniques used in this work and the 
possibility of improving on them. One is whether it is possible to extend the present 
result, or rather the general theorem used, to the case of non-constant equilibrium 
states, this is of vital importance if we want to consider non-flat backgrounds or 
even self-gravitating fluids. We think this is probable the case if we further assume 
that the theory has, at equilibrium, a conserved energy-that is a positive definite 
bilinear form in the tangent space of equilibrium states-as is usually the case for 
theories coming from a Hamiltonian formalism. The second possible extension is to- 
wards allowing non-compact Cauchy slices. There is another technique, introduced 
by Matsumura ||10|, which allows to study global existence and stability for some 



particular cases that ranges from hyperbolic heat conduction to relativistic super- 
fluids 0]. This technique allows to treat the case of non-compact Cauchy slices, 
but can not be applied to the general systems considered in this work. Thus, it 
would be important to extend the theorem given in the appendix B to the case of 
non-compact Cauchy slices. The third extension is in the direction of boundary 
values problems. It is clear that one would like to use this theory to describe situa- 
tions where the fluid is in a finite region of space, in that case the equations cease 
to be hyperbolic outside the region occupied by the fluid, and so a boundary value 
formulation is needed. Since Navier-Stokes fluids behave in a much more amenable 
way with regards to boundary conditions than perfect fluids, one would expect that 
these dissipative fluids will have that property too, making this an interesting, and 
perhaps tractable, problem. 
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APPENDIX A: EQUIVALENCE OF ASSUMPTIONS 

In this Appendix we prove the equivalence between the requirement (3) of Section 
2 and a condition imposed in literature 0, ^. We do this, for simplicity, in the case 
of divergence type fluid theories, and we assume that the reader is familiar with 
[0]. The condition under consideration allows to characterize the strict equilibrium 
states in such a way that they are precisely the same set of equilibrium states found 
for the standard Eckart theory. 

The requirement 3 of Section 2 is the following: The map J-'k : S^ -^ S* defined 
by J^xii'^) = KaM^AB i'^, is injective for all space-time vectors K"^ ^ 0, where 
S* denotes the dual of the space of fluid states and S^ the subspace of momentary 
equilibrium fluid states. 

The equivalence between this condition and the one assumed in [^ follows from 
the following argument. The map Tk is injective so, KaM^AB V'^ = ^ ip^ = 
where ip^ = {ip , ijj"' , 0) . Due to the definition of indices A and B, the system of 
equations above represents a scalar equation, a vector equation, and a symmetric 
two indices tensor equation. 

First, consider the scalar equation, the contraction of the vector equation with 
C, and the contraction of the two indices tensor equation with ("'(''. All this consti- 
tute a linear algebraic system of three scalar equations for variables Kaip""-, KaC°"ip, 
and 2KaC°'C^ipb- The injectivity implies that the only solution for these three vari- 
ables is zero and so, the determinant of the coefficient matrix is different from zero. 
Conversely, if the determinant of the coefficient matrix is different from zero, then 
we conclude that ■?/' = and ■j/'" = and then the map Tk is injective. By direct 
inspection it can be checked that these coefficients are the same found in equations 
(41)-(43) in 0. 

Second, consider the vector equation and the contraction of the two indices tensor 
equation with C^. This constitute a linear algebraic system of two vector equations 
for the variables mentioned in paragraph above, and 2C^'^K(^a'4'b) and -ft'fe'?/'. Because of 
injectivity, the only solution for all these variables is zero and so, the determinant of 
the coefficient matrix is different from zero. It can be checked, by direct inspection, 
that these coefficients are the same found in equations (45)-(46) in 0. 

Finally, consider the two indices tensor equation. It constitutes a linear algebraic 
two indices equation for the variables mentioned in the previous paragraph and 



K(^a'ipb)- Because of the injectivity, the only solution for all these variables is zero 
and so, the determinant of the coefficient matrix is different from zero. By direct 
inspection it can be checked that these coefficients are the same found in equations 
(48) in §. 

APPENDIX B: GENERAL STABILITY THEOREM 

Consider the Cauchy problem for a first order system of partial differential equations, 

^^ = {Ai + eA{iv,e))^+{Bo + eB,iv,e))v (Bl) 

v{t = 0) = fix^) 

where f : R x T'^ — > R", Ai{v,e) and Bi{v,e) are smooth (C°°) functions of their 
arguments, and f{x) : T'* -^ R" is also smooth. Let P denote the projector in the 
kernel of Bq. For the solution v of ( pT|) we define v^^^ = Pv{0, t) and w = v — v^^\ 



The stability theorem proved in M states|12 



Theorem 2 Suppose that the matrices Aq, Bq, and A\ are hermitian, and the sys- 
tem ^Bl\ j satisfy the '^relaxed stability eigenvalue condition", ie., the following con- 
ditions hold. 

(i) There is a constants > such that the eigenvalues X{k) of the symbol iAlkj + 
Bq satisfy Re{A} < —6 for all k & Q, k ^ 0. 

(a) The eigenvalues of Bq satisfy either Re{A(0)} < —6 or A(0) = 0. 

(Hi) ker Bq C ker Bi . 
Then, for < e < Eq with Eq small enough, the system I \B1\ ) is a contraction for w 
in a suitable norm, equivalent to a Sobolev norm H^ (p > s + 2), and v^^^ -^ const. 
when t ^ oo. 

The statement in the above theorem that the system is contraction means that there 
exists an if-norm, equivalent to the norm H^, such that 

j^\\w\\l<-{6 + 0{Emw\\j,, 

this implies that, for e small enough, there exists a global (in time) smooth solution 
of the Cauchy problem ( [Bl| ) such that there is an exponentially decaying bound 
for the Sobolev H^ norm of the w-part of its solution, and the v^^^ part goes to a 
constant when the time goes to infinity. 

The only difference between the fiuid equations (H) and (pl|) is the presence of N 
in front of the time derivative. This causes no difficulties, and the theorem above is 
applicable by a simple redefinition of the scalar product used. With this new scalar 
product, the eigenvalues problem for the fiuid equations becomes 



det (iN^^N^Qkj + Nq^Iq - \e) = 0. 



where E is the identity matrix. Then, the conditions (i), (ii), and (iii) of the theorem 
1^ can be verified as is done in the proof of theorem |l[ 
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